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X 1 Abstract 

We consider a quantum scalar field on an arbitrary gravitational background. 
We obtain the effective in-in equations for the gravitational fields using a 
covariant and non-local approximation for the effective action proposed by 
Vilkovisky and collaborators. From these equations, we compute the quan- 
tum corrections to the Newtonian potential. We find logarithmic corrections 
which we identify as the running of the gravitational constants. This run- 
ning coincides with the renormalization group prediction only for minimal 
and conformal coupling. 
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I. INTRODUCTION 



The Effective Action (EA) is an useful tool for analyzing the quantum corrections to the 
classical dynamics in quantum field theory. In particular, the effective equations derived from 
it should be the starting point to investigate many interesting problems like the influence of 
quantum matter fields on the behaviour of gravitational fields, both in cosmology and black 
hole physics. 

The EA is a very complicated object, even in the one-loop approximation, and it is 
necessary to develope approximation techniques in order to evaluate it. A widely used 
approximation is the Schwinger DeWitt expansion (SDW) 0, which consists basically in 
an expansion in derivatives of the background fields. This expansion is useful in situations 
where the background fields are slowly varying with respect to the mass of the quantum 
fluctuations and for the analysis of the renormalizability of the theory. However, many 
interesting physical effects are lost in this approximation. Alternatively, one can consider 
a situation where the background fields are weak but rapidly varying. In this case, it 
is possible to expand the EA in powers of the curvatures of the background fields. The 
resulting expansion has been recently investigated by Vilkovisky and collaborators 0- [|J, 
and it is in general a non-local object. 

On the other hand, there is a simple an intuitive way of taking into account, at least 
partially, the quantum effects. In quantum field theory, parameters like masses and coupling 
constants are not constants but scale-dependent quantities. This is due to vacuum polar- 
ization effects and the scale dependence is dictated by the renormalization group equations. 
One can use this fact to construct a 'Wilsonian' effective action ||, which is basically the 
classical action in which the parameters have been replaced by their running counterparts. 
An argument of this type has been recently proposed to explain the dark matter problem 
||: due to quantum effects the Newtonian potential should be modified according to 



where G(n) is the solution to the renormalization group equations in a renormalizable theory 
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(1) 



r 



of gravity with R 2 - terms in the Lagrangian |7j. As the theory is asymptotically free, G(r) 
is an increasing function of r, and this may explain at least part of the 'missing' mass. The 
running of G may also induce interesting cosmological and astrophysical effects || . 

It is the aim of this work to analyze the relationship between the non-local approach to 
the EA proposed by Vilkovisky et al, the renormalization group, and the 'Wilsonian' effective 
action. In Ref. ||, it was shown that the existence of non-local terms in the effective action 
is linked to the short distance behaviour of the theory and to the renormalization group. 
The analysis was done in a non-covariant weak- field approximation, at the level of the in- out 
effective action. Here we will extend that analysis: we will be using a covariant effective 
action, we will work at the level of the in-in (see below) semiclassical equations of motion and 
we will see explicitly the running behaviour of the gravitational constants in the Newtonian 
potential. For simplicity, we will consider a toy model in which we quantize a scalar field on 
a classical gravitational field. We will not include the graviton loop in our calculations. 

We would like to stress that our interest here is not to look for measurable corrections 
to the Newtonian potential. Indeed, we know a priori that these corrections are extremely 
small in the toy model considered. What we are going to discuss is how to derive from 'first 
principles' (i.e. the EA) the scale dependence of G in the gravitational potential. 

The paper is organized as follows. As a warm up, in the next Section we will obtain the 
running behaviour of the electric charge in QED starting from the non-local EA. In Section 
III we will consider a free quantum scalar field of mass m on an arbitrary gravitational 
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background. We will obtain the non-local version of the effective action in powers of — 
and the non-local effective equations. In Section IV we will compute the quantum correc- 
tions to the Newtonian potential and compare the results with the ones obtained from the 
Wilsonian' approach. Section V contains the conclusions of our work. 



3 



II. QUANTUM ELECTRODYNAMICS 



Due to the vacuum polarization, the electrostatic interaction potential between point 
charges is modified to 

e 2 (r) e 
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In the short distance limit (mr <C 1) we have 
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where r is defined by — In mr = 27 + | . 

On the other hand, the solution to the renormalization group equation gives the following 
running for the electric charge: 
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As can be easily seen from Eq. (|3]), in the short distance limit the electrostatic interaction 
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potential is just the usual ^ in which the electric charge has been replaced by its running 
counterpart Eq. (|]), with the additional rule that the mass scale fi is replaced by the inverse 
of the distance r and the mass scale reference [i is set at r~ l . Therefore the 'Wilsonian' 
argument gives the correct answer for the electrostatic potential in the short distance limit. 

We will derive here these old and well-known results using the EA formalism, since this 
exercise will be a useful guide to the more complex calculation presented in Sections III and 
IV. For simplicity we will quantize only the fermion field ( with no classical component ), 
keeping the electromagnetic field as a classical background. 

The classical action for QED in Euclidean space is given by 



S c = \f^ v + $(0 + ie4 + im)ifj 



(5) 



In the weak field approximation, the one-loop effective action obtained after integrating out 
the fermions is given by 
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s,„ = s c + ^j d >MA» (x )n Ax , xm + o<,A<) 



(6) 



where IT^ is the usual vacuum polarization tensor for QED. The renormalized effective 
action is 
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where 



F(n) 
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The modified Maxwell equations that derive from Eq. (0) are 
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where we included a classical source J. 
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The form factor F(D) admits the following integral representation in terms of the massive 
Euclidean propagator (M 2 — D) -1 : 



F(n) = i/ 1 A(1-^ 
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+ r 2 + (r=h 

Therefore, we can regard F(n) as a two point function whose action on a test function f(x 
is given by 



F(n)f( x ) = / d 4 x'F(n)(x,x')f(x') 



(11) 

All these equations are valid in Euclidean space. To get the Minkowski version of them, 
one should replace the Euclidean propagator by the Feynman one. However, the equations 
thus obtained are neither real nor causal since the effective action gives in-out matrix ele- 
ments instead of expectation values, making the interpretation of the equations awkward. 
Alternatively, one can use the Close Time Path (CTP) JO]] formalism to construct an in- in 
effective action that produces real and causal field equations for in-in expectation values 
0- 



The CTP formalism involves a doubling of the degrees of freedom and a generalization 
of the Feynman rules that includes both the Feynman and Dyson propagators, as well as the 
two-point Wightman functions (these functions carry the information about the quantum 
state of the system). However, if one is interested in the in — in effective equations for the 
standard m-vacuum state, this complication can be avoided. Indeed, in this situation it can 
be shown that the in — in version of the equations is obtained by replacing the Euclidean 
propagator by the retarded one |§ in the integral representation of the form factor Eq.[To|. 
Alternatively, the in — in form factor can be obtained in the in — out formalism by taking 



twice the real and causal part of the in — out form factor ||13|| . We will denote the in — in 
form factor thus obtained by F in (n). 

In particular, if the test function is time independent 

d 3 k 



F m (D)/(x) = F(V 2 )/(x) = J dV/^e tk -( x - x ')F(-k 2 )/(x') 



(12) 



because the time integral of the retarded propagator coincides with the Green function of 
the Laplacian. 

In the short- distance limit m 2 <C — □ the Euclidean field equation reduces to 
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One often encounters the distribution G(— -5 
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ln(— —), which will play a central role 



in what follows. The action of the in-in counterpart of G{ — on a test function f(x) is 



given by (see Refs (T^j , |L4| and also Hl5| ) 



G in (~)f(x) = '- I d'x'Oix - x°')5'((x - x'f)f{x') 
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where u and v are respectively the standard retarded and advanced coordinates with origin 
at the point x. When the test function is time independent, Eq. ( Ill ) reduces to 



G in (-£)/(x) = G(-^)/(x) 



(15) 



as expected from Eq. fll2). 

We are now ready to compute the modifications to the electrostatic potential. Taking as 
a classical source a static point charge, the modified Gauss law reads 
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The solution for the electric field is spherically symmetric E = E(r)f and we shall find it 
perturbatively in powers of e 2 



E = E (0) + E (1) 

The leading order term is the classical contribution 

V-E (0) = e5 3 (x) E <0 \r) 

and the first quantum correction is given by 
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Therefore we have to evaluate the action of G{—^r) on the delta function. Using the Eqs. 



(|P2|) and (JT5J) we readily obtain 
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where the last term gives a /i-dependent correction to the classical solution that will be 
absorbed into the classical source. The quantum correction is 
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where r is an arbitrary reference radius. Integrating Eqs. (|18|) and (^l|) and multiplying by 
the charge e we get the electrostatic interaction potential 
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From this equation the running behaviour of the electric charge follows, and it coincides 
with that of Eq.©. 



III. NON-LOCAL EFFECTIVE EQUATIONS FOR THE GRAVITATIONAL 

FIELD 



Let us now consider a quantum scalar field on a gravitational background. The Euclidean 
action for the theory is 



S — Sarav ~\~ S< 



grav > ^matter 



(23) 



where 
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+ m 2 (f) 2 + ^Rcj) 2 } 



(25) 



We have included terms quadratic in the curvature since in any case they will appear in the 
renormalization procedure. As we will use ^-function regularization, the constants G,A,a 
and j3 are finite (and dependent on a mass scale /z). 

The effective action for the classical gravitational field can be obtained by integrating 
out the quantum scalar field. Formally the result is 



-□ + m 2 + £R 
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def 



Sgrav r 



(26) 



where \x is an arbitrary mass scale. 

The evaluation of the above determinant in a general background is a very complicated 
task. Let us denote by 1Z either the Riemann tensor or any of its contractions with the 
metric. When the gravitational field is slowly varying, i.e. when V n TZ m <C m n+2m , one can 



use the SDW technique to get |L6 

1 
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where we have omitted /i-independent terms which redefine the classical constants. The 
functions cij (x) are the coincidence limit of the SDW coefficients, given by 

ao(x) = 1 

ai(x) = {±-t)R (28) 

a 2 (x) = —RuupaR^ - —RuuR^ - -(- - 0°R+ - (- ~ R 2 
' 180 M p 180 M 6 V 5 W 2 V6 V 



a n ( x ) = V 2n - 2 K + KSJ 2n - A K +■■■ + VVTZ^ 1 + IT (29) 

The last line shows schematically the coincident limit of a n (x). 

From the SDW expansion it is easy to derive the scaling of the gravitational constants. 
The effective action should not depend on the scale fi. As a consequence, taking //-derivatives 
in Eq. ( p6[) we find 

dG G 2 m 2 ( 1\ 
A*-T7 = — — U-« (30) 
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L L ~T~. = -7^Z-2 (32) 
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which is the usual running for the gravitational constants. We can use Eq. (|30|) to construct 
a 'Wilsonian' gravitational potential. The scaling of G is given by 

G^^Gofl + ^e-^ln^ (34) 

V 7T 6 fl J 

so the Wilsonian potential is 

In the next Section we will see if it is possible to derive this potential from the EA. 

The SDW expansion is not useful for the analysis of the short distance behaviour of 
the theory. As we have seen in the previous Section, one should consider weak but rapidly 



varying background fields. Assuming that W1Z 3> 7Z 2 , one may try to sum up in Eq. 
([27j) all the terms which contain a given power of the curvature. This rather complicated 
calculation has been performed by Avramidy in Ref. |3J]. See also Refs. || for the massless 
case. The result, up to second order in the curvature, is 
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Note the similarities between these form factors and the corresponding F(D) that appears 
in QED (Eq. (§)). 

From Eqs. ( p7[ ) and (|38|), one can derive the effective gravitational field equations. As 
we are neglecting 0(TZ 3 ) terms in the effective action, it makes no sense to retain 0(TZ 2 ) 
terms in the equations of motion. Therefore, when doing the variation of the action, it is not 
necessary to take into account the g^-dependence of the form factors. Moreover, it is possible 
to commute the covariant derivatives acting on a curvature, i.e., V ^S/ V TZ = V ' V V \JZ+0{1Z 2 ) . 
After a straightforward calculation we find 



1 mr tJ _ L. , , m 
8^ + l6^«-6 )( - 1+h V 

A m 4 . 3 , m 2 
+ -— -- + ln — 
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where 
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H$ = AV^V U R - Ag^UR + 0(R 2 ) 

H$ = 2W^7 U R - g^UR - 2UR pv + 0(R 2 



(40) 



and 
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F^HjV + F 2 (n)H 
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Up to here we made no assumptions about the mass m. In the large mass limit, m 2 lZ 3> 
WTZ the SDW expansion Eq. (p7|) is recovered (up to 0(TZ 3 )). However, as in QED, we 
are interested in the opposite limit. Let us assume that the typical scale of variation of the 
gravitational field is much smaller than m -1 , that is, m 2 TZ -C VV7£. In this situation, we 
can expand the functions Fi(n) and F 2 (D) in powers of — The result is 
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(It is possible to obtain exact expressions for F\ and F 2 in terms of elementary functions. 
However, we will not need these long expressions in what follows). Inserting the expansions 



Eq.(f42|) and Eq. (f43|) into the effective equations Eq.(^) we get 
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where we have set the scale \i so that the cosmological constant is zero and we included a 
classical source T^ s . 

As with the SDW expansion Eq. (|27|), one can easily derive //-dependence of the gravi- 
tational constants from these modified Einstein equations. Alternatively, as pointed out in 
Refs. |T7]] and ||, one should also see the running behaviour by performing the rescaling 



g^u —* s 2 9[j,v and looking at the large s limit. Since □ — > s 2 under this rescaling, the non- 
local terms proportional to In □ become relevant in this limit. From the terms independent 
of m in Eq. fl44j) we get 



i /.. i, 2 i 



a{ s) = a{ s=l)-— 2 ^--Y--yns (45) 

P(s) = j3(s = 1) In s (46) 

Hy ' HK ' 960tt 2 v ; 

It is worth noting that the scaling behaviour for a and j3 obtained using both methods Eqs. 

(HH) and (HID are identical. As far as the Newton constant is concerned, we can obtain 



its running behaviour only for £ = 0. In this particular case, the terms proportional to m 2 
in Eq. ( fHf ) have a logarithmic kernel that appears in the combination 

" >2 H-^)kHS>-2H$) (47) 



384vr 2 v m 2J U y 

Up to the order we are working (0(1Z 2 )), the basic tensors R^ v — ^Rg^ v , H^J and Hj 2 } are 
related by 

fl£> - 2H$ = 4D [R, v - ±Rg^ (48) 



Therefore, it is natural to express Eq. fl47|) using this relation, which thus leads to an s-scaling 
for G 

G(s) = G(s = l) (l - G{S = ^ )m2 In,) ,(^ = 0) (49) 



that is identical to the //-scaling Eq. (|34] ) for minimal coupling. For an arbitrary coupling 
we cannot get the s-scaling for G, since the logarithmic kernel does not appear in the simple 
combination Eq. ( P7Q . In this case the relation Eq. ( pE8"D makes the identification of the scaling 
behaviour ambiguous. 
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We shall see in the next Section how to obtain the running behaviour for G from the 
Newtonian potential. 



IV. QUANTUM CORRECTIONS TO THE NEWTONIAN POTENTIAL 

As in QED, the vacuum polarization effects contained in < T^ v > induce modifications 
to the Newtonian potential. We will now evaluate these corrections. To begin with, we 
must obtain the in — in effective equations. To this end, we should express the form factors 
as integrals of the massive Euclidean propagator and replace it by the retarded propagator 
(see Eq. (fLOf) ). However, when computing the Newtonian potential we will consider only 
time independent fields. Therefore, i^ n (D) = F(V 2 ) and the in — in equations are just the 
Euclidean equations with □ substituted by V 2 . 

In the static weak-field approximation we have 

R = ^V 2 h (50) 

where we assumed the Lorentz gauge conditions [h^ — \n^ iP h)^ v = 0. For a point particle 
with = 5®5® u M5 3 (x), the trace of the (linearized) Eq. © is 

1 .v 2 - 2(3a + 0) V 2 V 2 - ( ^- - ^ + 6£ 2 - 3(£ - \f H~X)\ V 2 V 2 



IQttG v 1 32tt 2 Vl80 3 v 6 7 v /i 
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+0(m 4 ) = -M5 3 (x) (51) 

For simplicity we shall compute only the trace h and not the complete h^, since this will 
simplify the calculations and will be enough for our purposes. In the limit a,/3 — > 0, —h is 
four times the Newtonian potential. 
We shall solve Eq.(O) perturbatively 

h = h {0) + (52) 
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The classical contribution h^ ' satisfies 

(V 2 - a- 2 V 2 V 2 )/i (0) = -167rGAf5 3 (x) (53) 
where a~ 2 = 327rG(3a + j3). The time independent and spherically symmetric solution is 



(54) 



The first quantum correction satisfies 
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We now find the solution to this equation. To begin with, we will consider the limit 
or — > oo, since in this approximation it is easy to find such a solution. In this limit the 
classical potential becomes 

/i (0) = AGM (- + 4vra- 2 (5 3 (x)") (57) 



Using the action of the kernel G{— ^-) on the delta function (Eq.(|2HD) we find 
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where the coefficients A, B and C depend on m, n and £. The solution to this equation is 

24G 2 Mm 2 1 , In £ 12G 2 M 
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(59) 



The first, second and third terms come from the sources r -3 , r~ 5 and r -7 . The dots denote 
a term proportional to the classical solution as well as corrections at the origin, which 
are proportional to <5 3 (x) and its derivatives. They all come from the sources proportional 
to A, B and C in Eq. (|5S|) . We have not included them because our quantum corrections are 
not accurate near the origin. Indeed, we have derived the modified Einstein equations under 
the assumptions W1Z ^> 1Z 2 and m 2 lZ -C VV-R. Both conditions are satisfied for the 
potential if GM C r <C m , so the origin r = is excluded. 



From Eq. (j59|) we see that there are two different types of terms in the quantum cor- 
rection. The term containing the logarithm comes from the non-local terms proportional 
to m 2 ln(— □) in Eq. (f44|) . It is qualitatively what we expected from 'Wilsonian' arguments. 
However, the coefficient MZiMzsf. (£ 2 _ _L) j s no t exactly the same as the one derived from 
the renormalization group equation ([34]), unless £ = or £ = |, i.e., minimal or conformal 
coupling. This is an important difference with respect to the QED calculation, and shows 
that the 'Wilsonian' arguments are not always quantitatively correct. Besides the running 
of G, we have found additional r~ 3 and r -5 corrections. 

There are no terms in that we could associate to a running of the constants a and (3. 
This is not surprising because such a running would imply terms of the form In ^5 3 (x), which 
are ill-defined. Moreover, we have already pointed out that our quantum corrections are not 
valid near the origin. Therefore, to see the running of these constants we shall evaluate the 
exact solution for h and then analyze the limit err — > (to this end it is necessary to consider 
only the case m 2 = 0) 

As a is proportional to \3ot + f3\~^lpj anck , this limit makes sense only for very large values 
of a and (3. Otherwise the limit would apply only for r smaller than the Planck length lpi an ck, 
where our semiclassical calculations are not valid. Experimentally |18| |o;|,|/3| < 10 19 , so 
r ^ ll lpianck- This is still an extremely small length. Therefore, what follows should be 
taken only as an exercise that shows that it is possible to extract the scale dependence of a 
and (3 from the potential h. 

The calculations for the exact solution to Eq.(|55D are presented in the Appendix. We 

15 



quote here the main results. In order to solve the linearized equation of motion we have to 
evaluate the action of the kernel G(—^-) on the Yukawa potential 

V72 -or 2 -or or -or 

G(--f)— = H- 2 )— - -Ei(-ar) - —Eiiar) (60) 

where Ei(x) is the exponential integral function. With the help of this formula we get the 
exact solution for h^ 1 ' (see the Appendix) and the limit ar — > can be taken. The solution 
reads 



7T 6 



3 

r In (err) + (7 — —)r 



0(er 6 ) (61) 



It is worth noting that the logarithmic term is exactly the one expected from the renor- 
malization group scaling of 3a + (3. Indeed, for small err the classical potential becomes, up 
to a constant, 

h(o) ~ -2MGa 2 r (62) 

Substituting in the above equation Ga 2 = (327r(3a + /3)) _1 by its running counterpart 
(Eqs.(|31~l-j32j)) with n = -, one finds 

llQ) 6G 2 Ma 4 L 2 , r 

hp) ~ ^__ 2 r ln— 63 

7r 6 r 

which coincides with the logarithmic term of result Eq. (|61|). 



V. CONCLUSIONS 

Let us summarize the new results contained in this work. We have obtained the in-in 
effective equations for an arbitrary gravitational field that include the backreaction produced 
by a quantum scalar field of mass m. The equations are non-local, covariant and valid under 
the assumption WTZ ^> 1Z 2 . In the limit m 2 ^> — □, the equations become local and 
reproduce the Schwinger DeWitt expansion. In the opposite limit, m 2 <C — □, the presence 
of non-local kernels of the form ln( — y) made it possible to read the scaling behaviour of 
the gravitational constants a and f3 under the rescaling of the metric. This scaling coincides 
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with the renormalization group predictions. This is also the case for the s-scaling of the 
Newton constant, but only for minimal coupling (this fact has been pointed out in Ref. ||). 

Using the in-in equations we computed the quantum corrections to the Newtonian po- 
tential. This is our main result. We have found two types of corrections: short range 
corrections that decay faster than - and corrections proportional to - In — , which we recog- 
nized as the scaling of the Newton constant. This scaling coincides with the renormalization 
group prediction only for minimal and conformal coupling. For other couplings, while the 
//-dependence of G is proportional to (£ — |), the scaling in the Newtonian potential is pro- 
portional to 6(£ 2 — t^j). Therefore, the 'Wilsonian' approach is strictly valid only for £ = 
and £ = |. 

One of the main motivations behind the present work was the remark made in Ref. 
about the possibility of explaining the dark matter problem through the scale dependence of 
G. In that paper, the running assumed was the one dictated by the renormalization group 
equations, in a theory of gravity containing i? 2 -terms. From our results we see that, in 
the toy model we have considered, the renormalization group behaviour is qualitatively but 
not quantitatively reproduced at the level of the Newtonian potential. However, at present 
we can not draw definite conclusions about the R 2 theory, since we have not included the 
graviton loop in our calculations. We hope to clarify this issue in the future. 

Finally, we would like to point out that the covariant effective equations we have found 
in Section III can also be used to analyze the effect of scaling in cosmological situations. For 
a Robertson Walker metric with scale factor a(t), we expect local terms of the form lna 2 (t) 
to be contained in the kernel ln(— t-). These local terms may have interesting cosmological 
and astrophysical consequences, like the generation of a primordial magnetic field during 
inflation, [plj Work in this direction is in progress. 



Note added: While we were writing this article we received a paper by Donoghue, |T9 
where the author calculates the quantum corrections to the Newtonian potential due to the 
graviton loop. His results are qualitatively the same as ours in the case m 2 = 0. This is to 
be expected, since the physical degrees of freedom of the graviton can be treated as massless 
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scalar fields. 
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APPENDIX A: 



In this Appendix we calculate the first quantum correction that solves the linearized 
equation of motion Eq. fl55p for the case m 2 = 0. The evaluation of the action of the kernel 
G(— ^r) on the Yukawa potential is accomplished using Eqs.([TJ) and ( |T5| ) 



d 3 k 



fk 2 ) e 



(27 



/i 2 r' 



p.2 g— or gcrr e _CTr 

In (— ) Ei(-ar) Ei(ar) 



(Al) 



ir r 

where Ei(x) is the exponential integral function. Taking this expression into account, the 
equation of motion reads 

5 



(v 2 -a- 2 V 2 V 2 ) /»«=X;/iW 



(A2) 



where 



/i(x 
/ 2 (x 
/s(x 
/ 4 (x 
/ 5 (x 



64G 2 Ma 2 
64G 2 Ma 2 
64G 2 Ma 4 
64G 2 Ma 4 
64G 2 M(7 4 



<5 3 (x) 



i=i 



19 5£ 3£ 2 3,. l. 2l 2 

- + — + -(£- -) In// 2 

1440 24 4 8^ 6 ; P 

— (£--) 2 

167T 6 jr 

19 5£ 3£ 2 3 L 2l a 

h — — H (£ - -) 2 ln- 

57607T 967T lQn 16n 6 /i 



32tt v 6 ; . 
32tt vs 6 ; . 



— Ei(-ar) 
r 

e -crr 

Ei(ar) 



(A3) 



Being a linear equation, we propose a time- independent, spherically symmetric solution 
of the form 

5 



(A4) 



i=l 



where each /&^(r) is the solution corresponding to the source /j(x). 

Let us denote by £?(x — x') the Green function of the operator V 2 — cr~ 2 V 2 V 2 (obviously 
G(x) is proportional to h^(x)). The solutions h\ 1 \'x) are then given by 



/ ij (1) (x) = /dV £(x-x')/;(x') 



(A5) 
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The source /i(x) is proportional to (5 3 (x). Therefore, hi is proportional to and can 
be absorbed into the classical parameters. 
For i = 2 we obtain 



ix z r o 



1 r<x> 
1 \2 



dt- 



\nt 



t(l + t 2 ) 

—G 2 Mo 2 {t - -f / dz\e z EA-z 
itr 6 Jo 



sin(ort) 

) ~ e- z E t {z)} 



(A6) 



where the last equality can be proved by taking r-derivatives on both sides and using prop- 
erties of Ei(z). Having now the exact first quantum correction, one can analyze the limit 



or — > 0. Using the series expansion for the exponential integral function ||20|| , the quantum 
correction reduces to 



7T 



rln(crr) + (7 — -j) rJ r 



0{o 6 



(A7) 



The other sources can be treated in a similar way. However, as they are all proportional 
to cr 4 (see Eq. (|A3|) ), the new solutions h\ % = 3,4,5 are of order cr 6 . Therefore, 



/i«(x) = /^(x) + O(o 



(AS) 
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